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ABSTRACT 

Two  methods,  the  generalized  moment  and  the  regression, 
both  based  on  the  limiting  distribution  of  the  multinomial,  are 
given  for  estimating  the  parameters  in  the  removal  trapping 
method  of  estimating  animal  and  insect  populations.  Some 
finite  sample  size  results  are  provided  indicating  the  speed  of 
convergence  to  the  limiting  distribution.  Numerical  examples 
are  also  discussed. 
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SIMPLIFIED  POINT  AND  INTERVAL  ESTIMATION 
TOR  REMOVAL  TRAPPING 

t 

Andrew  P.  Soms 


1.  Introduction 


A thorough  discussion  of  the  removal  trapping  method  of  estimating 
animal  and  insect  populations,  together  with  limitations,  is  given  in  [6], 
pp.  182-6.  It  is  pointed  out  in  [4]  that  this  method  is  particularly  suited 
for  insect  populations.  Briefly,  there  are  assumed  to  be  m organisms 
in  some  fixed  area,  k trapping  or  sweeping  periods,  k > 2,  and  each  organ- 
ism is  assumed  to  have  a constant  probability  p of  being  captured  in  any  of 
the  k periods,  independent  of  the  other  organisms  (The  organisms  are 
not  released  when  captured).  If  the  trapping  probability  is  p,  then,  as 
pointed  out  by  Moran  [5],  p.  308,  the  joint  density  of  the  n , 1 < i < k , 

the  number  of  organisms  trapped  in  each  of  the  periods,  is 
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where  0 < p < I . The  above  is  seen  to  be  a multinomial  distribution, 

with  k + 1 categories,  and  parameters  m and  p.  = p(l-p)  , 1 < i < k , 

k k 1 

and  p = (l-p)  = 1 - V p.  = 1 - V p(l-p)1'  . it  is  desired  to  estimate 

k+1  f i f 

m and  p and  give  asymptotically  exact  confidence  intervals.  In  [5]  a 
method  based  on  maximum  likelihood  is  proposed,  which  is  elaborated 
upon  by  Zippin  [7].  In  addition  to  replacing  m ! by  Stirling's  approxi- 
mation, the  effect  of  which  is  not  clear,  both  Moran  [5]  and  Zippin  [7] 
state  that  the  usual  regularity  conditions  for  the  joint  asymptotic  normality 
of  the  maximum  likelihood  estimators  are  not  satisfied  in  this  case  (in 
addition  to  other  assumptions,  it  is  assumed  that  the  parameters  being 
estimated,  m and  p,  remain  constant,  which  is  not  true  here,  since 
the  asymptotic  behavior  is  for  fixed  p as  m — *>),  and  then  they  pro- 
ceed in  the  hope  that  somehow  a justification  may  be  produced  without 
giving  it.  Further,  even  if  these  difficulties  are  neglected,  the  estimating 
equations  are  either  implicit,  requiring  iteration,  or  after  some  approxi- 
mations, require  charts.  Here  two  theoretically  justifiable  methods  are 
discussed,  both  based  on  the  limiting  distribution  of  the  multinomial, 
which  give  the  estimates  explicitly  as  functions  of  n^,  . . . , n^  . The 
first  is  a modified  method  of  moments  and  the  second  is  based  on  regression 
estimates. 
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2.  The  Modified  Moment  Estimates 


Since  En.  /En  = 1 - p = o,  2<i<k,  equating  expectations  to 

1 i-1  — — 

the  observed  values  gives  n /n.  ^ = I - p = q . Note  that  this  is  defined 

only  if  n , > n . To  minimize  this  effect  it  is  reasonable  to  take  the 
l-l  i 

geometric  mean  of  the  k - 1 estimates  to  obtain  the  estimate 


n? 

3 - i - ? = (-£ 

1 2 


nk  l/(k-l)  = \ l/(k-l) 
nk-l  nl 


and  this  is  the  estimate  to  be  considered  here  - note  that  it  fails  to  exist 


only  if  n > n , an  event  which  will  be  shown  to  have  limiting  probability 

K 1 

k k 

0 . Since  E(Z  n.)  = m(l-(l-p)  ),  the  moment  estimate  m of  m is 

k 

y n. 


m 


1 


l-(l-p)k 


Consider  now  the  problem  of  asymptotically  exact  confidence  inter- 
vals - a reasonable  assumption  is  that  p stays  constant  and  m -*■  so  . 
The  asymptotic  distributions  of  p and  m will  be  obtained  by  using  two 
results  - the  joint  asymptotic  normality  of  n = (n  , . • . , n ) and  a result 

I K 

given  in  Anderson  [1],  pp.  76-7.  It  is  well  known  that  as  (p  , . . . , p,  ) 


stays  constant  and  m -*  ao  , 


n -mp, 

( — T~,  1 <i  < k)  - N(0,R)  , 

(mp.q.)2 


(2.1) 


. w 


("*♦"  means  convergence  in  distribution),  where 
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R=fPij]’  Pii=1  andfor  l*i,  Ptj  = -(P.P./tQ.qj)) 


(see,  e.g.  Johnson  and  Kotz,  [3],  p.  284)  - recall  that  here  p.  = pd-p)1'1  , 
1 — 1 — k (also  for  notational  convenience  Pj  = P,  Oj  = 1 - p = q),  and 
hence  it  suffices  to  keep  p constant.  The  result  cited  in  Anderson  is: 

Let  f(x)  be  a function  with  first  and  second  derivatives  existing  in  a 
neighborhood  of  x - b,  b = (bf,  . . . , bfc)  a fixed  vector,  and  suppose 


I T W ~ 

V n (U(n)-b)  - N(0  , T)  . Then 


\'n(f(U(n))- f(b))  ^ N(0,$  ' TS>  ) , 

b b 


where  ® j,  - ( 9x  > • • • » 9x  )£  • The  asymptotic  distribution  of  1-p 

1 k 

will  now  be  obtained.  It  follows  from  (2.  1)  that 


(2.  2) 


N/m(n./(m(p.q.)?)  . (p./q.)^  1 < i < k)  - N(0,  R) 


(2.  3) 


where  R = [p..],  p..  = 1,  and  for  i * j,  p..  = -(p.^/fq.q.))?  . Take 

f(x ) to  be 


l/(k-l) 


xkv'pkqk 


(2.4) 


In  all  that  follows,  b = ((p/q^ (p^)’)  . Then,  using  (2.2), 


Vm((n  /n  )1/(k_1)  . f(b ))  = >/m((n  /n,)1/^"1)  . 


w 

(1-P))  - 


N(°’<CT>2<f  +k£fcr-^)  = N(0,,r2,  , 

pq  p 
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(|j-)£  = -( l/(k -1))  q(q/p)  2 , 

k -1  k -1  - 

(— )~  = (l/<k-l))  q((l- pq  )/(po  ))2  , 

k 


(2.  5) 


and  R<f>,  = o'  . Therefore  also 
b b p 


sfm  (p  - p)  - N(0,  <r  ) , 


where  p = 1 - (n  /n J1  -/(k'1)>  or,  equivalently,  p is  asymptotically 

K 1 


N(p,  cr  /m)  . Using  the  same  technique  on  m,  let 


k/(k-l) 


k 

g(x)  = (V  x,(p.q  ) ?)/ 

1 A 1 


1 - 


(2.6) 


Then 


8g 


( B^lb  = (piqi)2/(1-q  ) for  2 <i<k-l  , 


= ((pq)?  - (k/(k-l))qk(g/p)?)/(l-qk)  , 

oXj  b 


(li"  )£  = ((pkV' + (kAk-1))q(pkqk)2/p)/(1-qk) 

k 


w 


(2.7) 


Since  g(b ) = 1,  by  (2.  2) 


*s/m(m/m-l)^  N(0,  $^T$^)  , 


where  ® ‘ =(  )£  • In  order  to  evaluate 

1 k 


let  a'  = ((PjPj)2,  . . • , (Pkqk)  ) » 
b 
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k 1 1 

and  b'  = ( -(k/(k-l))(l-p)  (q/p) *,  0,  . . . , 0,  ( k/( k - 1) ) ( 1 - p) ( p q ) ?/p)  • Then 

$ ^T4>  = a'Ta  f 2a'Tb  f bTb,  and  after  some  algebra, 

b'Tb  = (k/(k-l)Mqk+1/p)(qk  + pqk_1  f 1)  , 
a'Tb  = 0 , 

k V ? 

a'Ta  = 1 - q - (1-q ) . 

Hence  = ( ( k /( k -1) ) ^ ( Qk +1/p) ( 1 + pQk “ 1 +q k ) + qk ( 1 - ak) )/( 1 - qk ) 2 = o-2  . 

u u m 

t A 2 A . ^ ^ 

So  v m (p-p)  — N(0,  o-  ) and  (m-m)/Vm  —N(0,cr  ).  Since  p—  p 

p m 

and  m/m  -*•  1 in  probability,  <r2  — <r2  , cr2  — cr2  in  probability  (o-2  and 

p p m m p 

(r*  are  obtained  from  cr2  and  o-2  by  replacing  p by  p,  q by  q ) and 
m pm 

A I A A i 

the  limiting  distribution  of  both  (p-p)/(o-A  /m2)  and  (m-m)/(cr~  m 2 ) is 

p m 

the  standardized  normal,  and  therefore  asymptotically  exact  (marginal) 

1 -a  confidence  intervals  for  p and  m are 


p ± z crA  /m  ? 

a P ’ 

ys  ^ 2 

m ± z a * m , 
a m 

th 

where  z is  the  upper  100a — percentile  of  the  standardized  normal.  A 
a 

Monte  Carlo  example  is  now  given  to  illustrate  the  asymptotic  theory. 

Three  cases  were  considered;  p = .4,  k = 3,  and  m = 100,  200,  and  400. 

In  each  case  1000  random  samples  (using  computer  generated  random  numbers) 
were  taken  and  for  each  sample  the  point  estimates  and  confidence  inter- 
vals calculated.  Two  coverage  relative  frequencies  were  computed  - the 
relative  frequency  Cl  of  the  estimate  lying  within  ± cr  and  ± 2c  of  the 
true  value,  and  the  relative  frequency  C2  of  the  estimated  confidence 
interval  covering  the  true  parameter  value  - the  latter  is,  of  course,  of 


L. 
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the  greatest  interest.  The  results  are  given  in  Table  1.  Here  a is  the 

A 

theoretical  standard  deviation  (st.  dev.  ) and  <j  is  the  standard  deviation 
estimated  from  a sample. 


1. 

Monte  Carlo  Results  for  1000  Samples 

a. 

P 

Sample 

Sample 

p 

m 

Mean 

cr 

st.  dev. 

Cl(±I(r) 

Cl(±2tr) 

C2(±la) 

C2(±2o-) 

. -1 

100 

. 397 

. 092 

• 090.. 

. 68 

. 96 

.70 

. 97 

. 4 

200 

. 399 

. 065 

. 68 

. 96 

. 68 

. 96 

. 4 

400 

. 401 

. 046 

. 04  5 

.70 

. 96 

.70 

. 96 

m 

. 4 

100 

105.  0 

1 3.  7 

22.  1 

. 69 

. 92 

. 77 

. 92 

. 4 

200 

204.  1 

19.  4 

22.  3 

. 69 

. 93 

.73 

. 93 

. 4 

400 

403.  7 

27.  5 

28.  8 

. 69 

. 94 

. 69 

. 94 

An  alternative  procedure  to  using  a-  in  the  confidence  intervals  is 
to  correct  for  bias  in  the  point  estimate  (e.g.,  for  m = 100,  5.0  is  sub- 
tracted from  the  estimate  and  this  is  considered  to  be  the  new  estimate) 
and  to  use  the  sample  standard  deviation  s (here  based  on  1000  random 
samples).  As  might  be  expected,  this  results  in  slightly  conservative 
intervals.  The  results  are  given  in  Table  2,  based  on  1000  simulations 
in  each  case.  In  practice  the  sample  estimates  would  be  used  as  the  true 
values  in  the  simulation. 
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r 


2.  Coverages  with  Bias  Substracted 


when 

Sample  Standard  Deviation  is 

U sed 

a. 

P 

p 

m 

Cl(±s) 

Cl(±2s) 

C2(±s) 

C2(±2s 

. 4 

100 

. 68 

. 96 

. 70 

. 97 

. 4 

200 

. 68 

. 96 

. 68 

. 96 

. 4 

400 

. 70 

. 96 

. 70 

. 96 

b. 

m 

. 4 

100 

. 61 

. 93 

. 89 

. 97 

. 4 

200 

. 6 5 

. 94 

. 73 

. 96 

. 4 

400 

. 67 

. 95 

. 69 

. 96 

Sometimes  before  the  trapping  experiment  is  begun,  a preliminary 
estimate  of  p is  available.  In  this  case  k,  the  number  of  trapping  periods, 
can  be  chosen  so  as  to  minimize  the  variance  of  p — because  of  the  im- 
portance of  p in  this  method  of  estimation  this  is  a reasonable  optimality 

criterion.  The  function  to  be  minimized,  apart  from  multiplicative  con- 

2 k-1 

stants,  is  f(k)  = l/(k-l)  (1+1/q  ) . The  minimum  is  obtained  by  setting 

the  derivative  equal  to  0,  the  resultant  equation  being 

qZ  + 1 = (-fn  q/2)z  , 

z = k - 1 . Table  3 gives  the  nearest  integer  to  the  exact  minimum  value, 
which  is  readily  obtained  by  iteration,  as  a function  of  selected  q . 
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3.  Optimum  k Values 


q_ 

95 

90 

85 

80 

70 

60 

50 


k_ 

44 

22 

15 

11 

7 

5 

4 


The  customary  statistic  used  to  test  the  adequacy  of  the  model  is 


Z = ^ (n.  - mp.)2 * 4 *,/(mp.),  where  p.  = p(l-p)i  . It  is  not  at  all  clear, 
in  this  or,  of  course,  the  maximum  likelihood  case,  that  Z has  an 


asymptotic  (p  fixed,  m -*■  =o)  ^ distribution  with  k -2  degrees  of 

freedom  (d.  f.  ),  since  the  usual  regularity  conditions  (see  [2],  pp.  500-1, 

506)  are  not  satisfied.  The  empirical  approach  given  here  consists  of 

using  (2.  2)  to  obtain  the  expected  value  of  the  limiting  distribution  of  Z 
2 

and  then  to  fit  a \ distribution  (as  is  done  with  good  results  in  fitting 

2 

the  distribution  of  sums  of  x random  variables)  by  estimating  the  d.  f. 
using  the  parameter  estimates.  The  observed  value  of  Z is  then  compared 


th  2. 

to  the  upper  100a — percentile  of  the  fitted  x (using  interpolation  on  the 


d.  f.  , since  in  general  the  fitted  d.  f.  will  not  be  integral).  Specifically, 
consider 


4 — — — i - 1 

f (x)  = x.  (p.q  ) 7 - g(x  )(1  - f ( x ) ) (f(x  ))  , 


l l l 


J 
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where  f(x)  end  g(x  ) ire  given  by  (2.  i)  and  (2.6),  respectively. 
Note  that  f.(b)  0 and  thus  from  (2.2), 


\'mf.  (n  ''(m(p,q  )'),  . . . , nk  /(m(Pkqk)  *))  - si  mf.  iu  ) 


— aaa[.|  w 2 

= vm((n  - mpq  ) m)  -*  N(0,  a . ) , 


(2.  8) 


&f,(x)  9f.(x) 

l l 


where  <r.  = <&  ;~  T$  ~ = ( — - , , 

\ lb  lb  lb  dx 


9x  'b 
k 


Using  (2.  5)  and  (2.7)  and  the  chain  rule,  an  expression  for  <r.  can  be 


obtained,  e.g., 


9f.(x .)  i n ~ . , 

<-fcr->E  - - hT’S  •*.  “‘S'"*  • 

i i 


and  similarly  for  the  other  derivatives.  Even  though  explicit  expressions 

2 

do  not  appear  practical,  o\  is  easily  evaluated  by  means  of  a short 

i 

computer  program.  Since  (m/m)?  -►  1 and  p.  -*  p in  probability,  it 
follows  from  (2.  8)  that 

a ^^i -1  , ~ A 2 . 

(n.  - mpq  )/(mp.)  - N(0,  <r  . /p. ) , 


vm  /N  | _ 1 2 /v/v 

and  therefore  the  asymptotic  mean  of  Z = ' (n,  - mpq  ) /(mp.)  is 

k 2 f 1 1 

p = Vo-/  p,  . Replacing  the  parameters  p.  and  q.  by  their  estimates 

j i i ii 

a /v 

p.,  q in  p gives  the  estimated  d.  f.  of  the  distribution  of  Z,  and  using 
these  d.  f.  a cut-off  point  for  the  adequacy  of  fit  test  can  be  obtained 
from  tables. 
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3.  Regression  Estimates 


The  method  of  moments  given  above  depends  on  the  "extreme" 

values  n and  n,  and  therefore  in  cases  where  the  trapping  periods 
1 k 

are  too  short  and  hence  the  data  show  considerable  fluctuation,  it  may  be 
preferable  to  use  an  estimate  of  p depending  on  all  the  data.  As  sug- 
gested in  [4],  a simple  check  is  to  plot  n^  against  i on  semi-log  paper 
and  see  whether  the  plot  approximates  a straight  line.  If  yes,  the  method 
of  moments  with  its  attendant  computational  simplicity  can  be  used,  and 
if  no,  then  the  regression  method  discussed  below  may  be  used. 

It  is  pointed  out  in  [5]  that  log  En^  = log  +u  lie  on  the  straight  line 

log  p.  = i log(l-p)  - log(l-p)  + log  p + log  m , (3.  1) 

but  this  method  then  is  dismissed  in  [5]  by  saying  that  the  usual  assump- 
tions of  regression  theory  are  not  satisfied.  Here  a different  approach  is 
taken  - namely,  the  point  estimator  of  p suggested  by  regression  theory 
is  used  but  then,  in  place  of  the  usual  regression  theory,  the  limiting 
distribution  is  obtained  from  (2.  2).  The  regression  equation  suggested 
by  ( 3.  1)  is 


log  n^  = ip  + a + e^,  1 < i < k , 

where  p = log(l-p)  (to  any  base),  a a constant,  and  ei  the  error  term 

A 

which  will  be  of  no  interest  here.  The  least  squares  estimate  p of  p is 


k 

p = V (log  n )(i 


k+1 

2 


k 


■>/£(* 

1 


k+1  2 
2 ' 
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i 


Since 


V 4 

- J 

1 


i = (k)(k+l)(2k  + l)/6  , 


p = V (log  n.)(i  - —j— )/(k(k*~  -1)/12) 
1 1 


Then  the  corresponding  estimator  1 - p of  1 - p is 


k c. 
1-P  = TT  n.  1 

i 7 


(3.2) 


where  c.  = (i  - (k+l)/2)/(k(k  -1)/12)  . In  order  to  show  that  ( 3.  2)  is 
consistent  for  p and  to  obtain  its  asymptotic  variance,  it  is  just  as 
easy  to  consider  general  estimates  of  1-p  of  the  form  (3.  2)  with  the  c. 
arbitrary  and  to  determine  the  conditions  on  c.  needed  for  consistency. 
Let 


h(x ) = TT (x  (p  q ) ?)  1 . 

1 11 


(3.  3) 


i-1 


Then,  using  (2.2)  and  p.  = p(l-p)  , 


lc  C*  ^ lx-  f~' 

\Tm  (TT  n '/m  - TTp.S 

1 1 


k k k-1 

k r H c-  lc'  Sic-i 

— i 1 1 f 1 V 1+1 

\frn  ( | | n.  /m  - p (1-p) 

1 


) • 


Therefore  a sufficient  condition  for  consistency  is  ^ c,  = 0 and 

k-1  1 1 

^ ic  , = 1,  and  in  this  case 
^ i+1  ’ 
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where  cr^  is  determined  in  the  usual  way  from  (3.  3)  using  (2.2).  It  is 
P 

noted  that  for  c.  = (i  - (k+l)/2)/(k(k? -1)/12)  these  two  conditions  are 

1 k 2 
satisfied,  since  clearly  V c = 0 and,  letting  c = k(k  -1)/12  , 

1 


k 1 k-1  k-1  k-1 

v 1C.  , = v 1(1+1  - ti!)/c  = (£  i2  - (k-1)/2  E ‘l/c 

f 1+1  t 1 1 


= ((k-l)(k)(2k-l)/6  - k(k  -l)/4)/c  = 1 . 


So,  for  this  choice  of  c , cr  = , 

1 p b b 


•i,  ■ ,xi(piqi,t)  ‘ - 


i 


(3.4) 


Therefore,  after  some  algebra,  = q2(T,  Ci(qi//pi)  ‘ 2 ZZ  CiCj)  ‘ 

A satisfactory  estimator  of  m is  obtained  by  the  same  argument  as 
for  the  moment  estimator,  namely, 

m = V n /(I  - (1-P)k)  , 

1 1 

where  1-p  is  now  given  by  (3.  2).  As  before,  let 

k k cj  ]< 

g(x ) = Yj  A1  - (TT (x.  n^pToT  ) ) ) • 
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Then,  by  (Z.  2), 


-1)2  N(0,a2m)  , 


Z 9g  9g  , , 

where  <r  = <*  1 T<i> * 1 = ( — , . . . , 5-  )£  and 

b b b 1 k 


I Ik 

(P.qj2  + kc.(q./p.)?  (1-p) 
ag  1 1 1 1 1 

3x.  'b  , „ ,k 

1 1 - (1-P) 


3.  5) 


^ k ^2  Zk  k ^ q. 

Therefore,  after  some  algebra,  <r  = ^ ~ ( V c,  — — - 

’ y ’ m . k k 2 ^ 1 p, 

1-q  (1-q  ) 1 1 

Z V c c ) . 

1 < j / 1 /( k - 1 ) 

It  should  be  noted  that  the  moment  estimator  of  1-p,  (n  /n  ) ‘ , 

is  a special  use  of  (3.  2)  with  c = -l/(k-l),  c.  = 0,  Z < i < k-1,  and 

1 k-1  1 

c^  = l/(k-l),  since  c^  + c^  = 0 and  V ic.+^  = (k-1)/ (k-1)  = 1 . 

The  discussion  of  lack  of  fit  is  completely  analogous  to  the  method 
of  moments,  except  that  the  partials  are  different  but  again  obtained  by 

the  chain  rule  using  (3.4)  and  (3.  5)  and  evaluated  by  a short  computer 

— I ~ ~ — i - 1 

program.  Specifically,  let  f.(x)  = x,(p.q.)2  - g(x  )(l-h(x  ))(h(x  )) 


1 1 1 


Then 


9f. 


(a“)u  = P^1"1  - (i-l)PQ* 


9x.  b 
J 


9x,  'b 
J 


i-2,9h(x)  i-1  9h(x) 

9x.  'b  4 1 9x.  'b  * 

J ] 


for  i * j,  and  if  i = j,  (p.q.)?  is  added  to  the  above  expression  in 
which  j has  been  replaced  by  i . 
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4.  Numerical  Examples 


The  two  methods  of  estimation  discussed  above  are  applied  to  the 
data  in  [4]  and  [5].  First,  consider  the  rat  data  in  [5]  consisting  of 
k = 18,  n = (49,  32,  31,  34,  16,  33,  22,  27, 17,  19, 18, 16,  18,  12, 14,  12, 17, 7)  . 

If  a semi-log  plot  of  n against  i is  made,  it  can  be  seen  that  there  is 
a large  variability  in  the  data,  suggesting  that  the  intervals  are  too  short. 

In  addition,  n appears  large  and  n small  compared  to  the  line  sug- 
gested  by  the  other  points,  and  thus  the  moment  method  should  under- 
estimate m,  which  Table  4 shows  (This  is  also  reflected  in  the  large 
calculated  Z value).  If  however  the  data  are  grouped  into  9 trapping 

periods  with  n . = n , + n i = 1,  2,  . . . , 9,  the  semi-log  plot  is 
gi  2 i - 1 2i 

much  smoother  and  the  moment  method  gives  comparable  point  estimates  to 
the  regression  method,  even  though  the  variances  are  somewhat  larger,  as 
is  seen  in  Table  4.  The  maximum  likelihood  estimates  and  their  estimated 

A A 

standard  deviations  (s.  d. ' s)  given  in  [5]  are  m = 520,  32.  9 and  p = . 07  56  , 
. 00933  (the  estimated  standard  deviations  for  the  estimates  discussed  here 

. A i A A 4 A 

are  understood  to  be  j*  /m  2 for  p and  a * m d for  m).  Note  that  in 

p m 

terms  of  p for  the  original  data,  the  probability  of  capture  p'  for  the 
grouped  data  is  p'  = p + p(l-p)  . 


-1.  Moment  and  Regression  Estimates  for  the  Rat  Data 


Original  data  Grouped  data 


Parameter 


Moments 

Regression 

Moments 

Regression 

p 

. 108 

. 07  37 

. 14  1 

. 139 

s.  d. 

. 02  12 

. 0100 

. 0260 

. 0186 

m 

4 52 

527 

529 

533 

s.  d. 

29.  4 

37.  0 

51.  0 

39.  0 

Test  of  fit 

Z 

31.  9 

17. 4 

2.  00 

2.  08 

d.  f. 

23.  3 

16.  1 

11.  2 

7.  11 

It  is  pointed  out  in  [4]  that  removal  trapping  is  sometimes  the  only 
feasible  method  of  estimating  insect  populations.  The  example  discussed 
there  consists  of  the  number  of  maccolaspis  flavida,  n.,  caught  in  k = 10 
sweeping  periods  (the  data  in  [4]  were  actually  grouped),  with  n = 

(72,  63,  44,  32,  31,  23,  17,  18,1  1,  13).  Using  a somewhat  involved  graphical 
method,  which  does  not  yield  any  interval  estimates,  Menhinick  obtains 
p = .212  and  m = 359.  The  complete  results  for  this  example,  for  both 
the  moment  and  regression  method,  are  given  in  Table  5. 
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5. 


Moment  and  Regression  Estimates  and  Fitted  Values 
for  the  Insect  Data 


r 


Parameter 


Moments 

Regression 

p 

. 17  3 

. 188 

s.  d. 

. 0289 

.0198 

m 

381 

37  0 

s.  d. 

24.  8 

14.  8 

Test  of  fit 

Z 

4.  59 

3.  28 

d.  f. 

12.  3 

8.  29 

A 

Data 

Predicted  Values  n. 

l 

A 

A 

A 

A 

72 

n. 

i 

n.  -n. 
i l 

n. 

i 

Vni 

66.  0 

-6.  0 

69.  5 

-2.  5 

63 

54.  5 

-8.  5 

56.  5 

-6.  5 

44 

45.  1 

1.  1 

45.  9 

1.  9 

32 

37.  3 

5.  3 

37.  3 

5.  3 

31 

30.  8 

-.  2 

30.  3 

- .7 

23 

25.  5 

2.  5 

24.  6 

1.  6 

17 

2 1.  1 

4.  1 

20.  0 

3.  0 

18 

17.4 

-.  6 

16.  2 

-1.7 

11 

14. 4 

3.  4 

13.  2 

2.  2 

13 

11.  9 

-1.  1 

10.  7 

-2.  3 
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5.  Concluding  Remarks 


The  purpose  of  this  paper  has  been  to  give  two  statistically  justifiable 
and  computationally  simple  methods,  the  moment  and  the  regression,  as  an 
alternative  to  the  maximum  likelihood  approach  which  suffers  from  two 
deficiencies:  the  standard  regularity  conditions  for  the  joint  asymptotic 
normality  of  the  maximum  likelihood  estimators  are  not  satisfied  and  the 
estimating  eauations  are  either  implicit  or  require  the  use  of  charts. 

Listings  of  short  computer  programs  that  calculate  either  the  moment 
or  regression  estimates,  standard  deviations,  and  the  adequacy  of  fit 
statistics  are  available  from  the  author.  The  moment  and  regression  estimates 
and  their  standard  deviations  are  also  readily  computed  by  hand. 

The  methods  discussed  here  should  also  be  useful  in  other  cases  where 
the  data  is  multinomial  and  the  standard  maximum  likelihood  regularity 


conditions  are  not  satisfied. 
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18  supplementary  notes 


1 19.  JEEV-WORClT7ctT»RnPAO'Tr«v*f»e*ldAenoTAA**Ty  end  Identity  by  block  number) 


Removal  trapping,  moment  estimator,  regression  estimator,  asymptotic 
confidence  intervals 


20.  ABSTR*^T  (Continue  on  reveree  elde  U noceeeery  and  Identify  by  block  number) 

Two  methods,  the  generalized  moment  and  the  regression,  both  based 
on  the  limiting  distribution  of  the  multinomial,  are  given  for  estimating  the 
parameters  in  the  removal  trapping  method  of  estimating  animal  and  insect 
populations.  Some  finite  sample  size  results  are  provided  indicating  the  speed 
of  convergence  to  the  limiting  distribution.  Numerical  examples  are  also 
discussed. 
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